Evidence for impurity-induced frustration in La2Cu04 
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Zero-field muon spin rotation and magnetization measurements were performed in 
La2Cui_3;M^04, for < a; < 0.12, where Cu^+ is replaced either by M=Zn^+ or by M=Mg^+ 
spinless impurity. It is shown that while the doping dependence of the sublattice magnetization 
{M{x)) is nearly the same for both compounds, the Neel temperature {Tn{x)) decreases unam- 
biguously more rapidly in the Zn-doped compound. This difference, not taken into account within 
a simple dilution model, is associated with the frustration induced by the Zn^"*" impurity onto the 
Cu^"*" antiferromagnetic lattice. In fact, from Tn(x) and M{x) the spin stiffness is derived and found 
to be reduced by Zn doping more significantly than expected within a dilution model. The effect of 
the structural modifications induced by doping on the exchange coupling is also discussed. 
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La2Cu04 has been the subject of an intense re- 
search activity since after the discovery of high temper- 
ature superconductivity (HTcSC).^ In fact, besides be- 
ing the parent of HTcSC it has been early realized that 
La2Cu04 is one of the best prototypes of two-dimensional 
S=l/2 Heisenberg antiferromagnet on a square lattice 
(2DQHAFSL)i^ Accordingly, an impressive amount of 
theoretical and experimental studies on the effects of 
charge and spin doping in La2Cu04 have foUowedi^ 
In this respect a significant interest has attracted Zn^'^ 
S=0 for Cu^+ S=l/2 substitution, which is known to 
be one of the most detrimental for HTcSC— 1^ It has 
been observed that in La2Cu04 the substitution of Cu^+ 
with a spinless impurity can still be basically described 
by the 2DQIIAFSL Hamiltonian, with a spin stiffness 
(ps) renormalized by the spin dilution.—^— Namely, for 
La2Cui_a;Znj;04 one has 



H = -jY,Po.{x)Si.p^{:x)Si 



(1) 



where J is the superexchange coupling among the near- 
est neighbour (n.n.) spins and Poo{x) is the proba- 
bility to find a spin at site i or j for a doping level 
X, which for x — 5- can be approximated to (1 — x). 
Then, the above Hamiltonian can be mapped onto a 
2DQHAFSL Hamiltonian, with 2'Kps[x) = 1.15J(1 - 
a;)^i^il This simple dilution model is able to describe 
to a reasonable extent the basic properties of La2Cu04 
doped with S=0 impurities and, in particular, the main 
features of the zero-temperature sublattice magnetiza- 
tion M{x) and of the Neel temperature Tn{x)^^— Nev- 
ertheless, more recent accurate studiesiiii^ have clearly 
pointed out that M{x) and Tm{x) experimental data for 
La2Cui_2:Zn2;04 do not agree with the most accurate de- 
scription of a 2DQHAFSL achievable through numerical 
simulations. It is concluded that, although the dilution 
model is a good starting point to analyze the properties 
of La2Cui_a;Zn2;04, a complete understanding of the ef- 
fect of a spinless impurity in La2Cu04 is still missing. 



Two years ago Liu et al..— have suggested that S=0 im- 
purities might induce magnetic frustration. In fact, if the 
energy of the lowest unoccupied impurity orbitals is close 
to the one of Cu^"*" 3dx2-y2 orbital, sizeable next n.n. su- 
perexchange couplings may arise and compete with the 
n.n. coupling. ^"^ Accordingly, a reduction of Ps{x) and of 
Tn (x) faster than expected according to a dilution model 
is envisaged. 
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FIG. 1: Typical temperature dependence of the macroscopic 
static uniform spin susceptibility (closed squares, right ver- 
tical scale) and of the magnetic volume fraction as derived 
from /iSR measurements (open squares, left vertical scale) in 
La2Cui_2,Mgj;04. 



In order to test if impurity-induced frustration is re- 
ally relevant we have performed a systematic compar- 
ison of AI{x) and Tjq{x), determined by muon spin 
rotation (/xSR) and magnetization measurements, in 
La2Cui_j;Zna;04 and in La2Cui_a;Mg2;04, for equal dop- 
ing levels in the range < x < 0.12. Since the energy of 
the unoccupied Mg^"*" orbitals is far away from the one 
of Cu^"*" 3d^2_y2 orbital, no competing superexchange 
coupling should be present for Mg-doping, while for Zn- 
doping sizeable next n.n. superexchange couplings may 
arise. It is found that La2Cui_a;Mga;04 can be described 
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FIG. 2: Time evolution of tlie zero-field muon asymmetry for 
the X — 0.05 Zn-doped (top) and Mg-doped (bottom) sam- 
ples, at T = 50 K. The solid lines are the best fits according 
to Eq. 2 in the text. 



within a dilution model provided that minor modifica- 
tions of J due to the structural deformations are consid- 
ered. On the other hand, in La2Cui-a;Zn:r04 a reduc- 
tion of ps (x) faster than expected for a dilution model is 
found. We argue that this decrease of ps (x) indicates that 
S = impurities induce frustration^^ It is also pointed 
out that for Zn doping the behaviour of M{x) cannot be 
simply described by a 2DQHAFSL Hamiltonian, since 
the geometry of the underlying magnetic lattice is pro- 
gressively changing from the square lattice one to the one 
of a diluted J1-J2 model"'^'* as the doping increases. 

Polycrystalline samples of La2Cui_iM2:04 (M=Mg or 
Zn, x=0'02, 0.05, 0.075, 0.1 and 0.12) have been ob- 
tained by standard solid state reaction of 99.99% pure 
CuO, La203 and MgO or ZnO. Stoichiometric quanti- 
ties of these oxides were thoroughly mixed, pressed into 
pellets and multiple sinterings (12 to 24 h), with inter- 
mediate grindings, were performed under air, in the tem- 
perature (T) range 900-1100°C. In order to avoid oxygen 
excess, a last thermal treatment of 12h was performed 
under Ar gas fiow at 800°C. X-ray diffraction showed 
that the samples were single phase. Rietveld analysis was 
then performed using the space group Bmab and the lat- 
tice parameters were obtained. The T-dependence of the 
static uniform spin susceptibility, determined with a su- 
perconducting quantum interference device magnetome- 
ter, is characterized by a peak which marks (Fig- 1)- 
The sharpness of the peak confirms the correct oxygen 
stoichiometryji^ 

Zero-field (ZF) /j,SR measurements were performed at 
ISIS pulsed muon source on the MUSR beam line. Below 
Tjv the ZF muon polarization was characterized by clear 



oscillations associated with the onset of the long-range 
magnetic order (Fig. 2) and by the correspondent decrease 
in the longitudinal component of the muon asymmetryfi^ 
It was noticed that, particularly for La2Cui_xMgx04, 
two muon precessional frequencies are evident (Fig. 2), 
indicating two different muon sites. Thus, the decay of 
the ZF muon asymmetry could be nicely fit tcJ^ 



A{t) = Aie-'''*cos{jf,B'^t + 0) + 
+A2e-'''''cosi-f^B!^t + (j)) + Ase-^* + Bck , 



(2) 



where 7^ is the muon gyromagnetic ratio, B^2 is the 
local field at the muon sites 1 or 2, which are charac- 
terized by different dipolar hyperfine couplings Di yield- 
ing i?^ — Di{S), with (S) proportional to the sublattice 
magnetization. 171^2 are the decay rates of the oscillating 
components, which increase with doping owing to the 
enhancement of the local field distribution with increas- 
ing disorder. Bck is a small constant term arising from 
muons stopping in the sample environment and A is the 
decay rate of the muon longitudinal polarization. While 
in La2Cui_a;Zn3;04 Ai :S> A2, in the Mg-doped system 
the amplitude of the two components is similar (Fig. 2). 
The different ratio A2 j A\ found for the two systems in- 
dicates that, in spite of the same ionic charges and of the 
nearly equal ionic radii of Zr?^ and Mg^^ , the lattice 
potential around the two types of impurity is different. 
From the drop in the longitudinal component A3 of the 
muon asymmetry below Tjv it is possible to derive the 
magnetic volume fraction Vm, since for a powder sample 
V,n = (3/2) [1 - (^3/(^1 +A2+ Aa))].!? One observes 
a sharp increase in Vm at Tn (Fig.l). In Fig. 3 the T- 
dependence of i?^, normalized to its value extrapolated 
to T — >■ i?i (0), is reported as a function of T/T/v. One 
notices that, no matter which is the doping level, the T- 
dependence of the local field at the muon is the same once 
the two aforementioned quantities are properly rescaled. 

In Fig.4 Tn{x) and M{x)/M{0), derived from the ra- 
tio B'^{x,T 0)/Bf(0,r 0), are reported for the 
two systems. First of all, one notices that Zn doping 
induces a more rapid decrease of Tjv than Mg substitu- 
tion. In both cases Tn{x) shows a linear decrease de- 
scribed by Tn{x)/Tn{0) = (1 — ax), with an initial slope 
ctzn = 3.52 ± 0.17 for Zn doping and umq = 2.7 ± 0.15 
for Mg doping. Given the energy separation between the 
first excited Mg^+ orbitals and Cu^+ id orbitals, the lat- 
ter a value should be taken as the one expected for a 
diluted 2DQHAFSLii. (a = 3.196), reduced by the inter- 
plane coupling J_Liii However, to properly analyze those 
results, one has first to consider if the lattice deforma- 
tion induced by doping may cause a significant varia- 
tion in the n.n. coupling J.l^ Namely if J in Eq. 1 
depends on doping. First of all we point out that the 
changes in the lattice parameters show a similar trend 
for both types of substitutions (Fig. 3). One notices also 
that the c axis changes less in the Zn doped compound, 
where Tm decreases faster. Moreover, even in the x=0.05 
samples where the difference in La2Cui_2;Zna;04 and 
La2Cui_2:Mg2:04 latticc parameters is less than 1/1000 
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FIG. 3: (Top) The local field at the muon site 1, rescaled by 
its value for T ^ 0, is reported as a function of T/Tn{x) in 
La2Cui_2,Zna;04. (Bottom) Doping dependence of the c-axis 
length in La2Cui_2,Mgi:04 and in La2Cui_a:Zna;04. 
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a significant difference in the Neel temperatures is ob- 
served. These latter observations show that the major 
differences in the Tn{x) of the two compounds should not 
be associated with the effect of the lattice distortions on 
J{x) but rather, as will be discussed later on, with the 
frustration induced by Zn-doping. 

One can estimate the variation of J{x) due to the lat- 
tice strains on the basis of the experimental data for 
Tm{x) and M{x) reported in Fig. 4. In fact, one can 
write^'i'i^ 



Tn{x) = J^Pooix) 



fMjxW- 
\M{0)) 



(3) 



Since Tn ^ J one has £,{Tn,x) ~ exp{2TTps{x)/TjM{x))? 
Now, assuming that for La2Cui_2:Mg2:04 a simple dilu- 
tion model works one can write 2ttp^s{x) = 1.15J(a;)(l — 
x)^ (see Eq.l and RefH). For a fixed x value one notices 
from Eq. |3]that an increase in J{x) causes a correspon- 
dent enhancement of T/v (x) . By taking the experimental 
values for Tn{x) and M{x) one derives an increase of 
J{x) by only 40 K for the x=0.12 sample^^ with re- 
spect to J(0) ~ 1580 K?2, This increase is close to the one 
that could be estimated from the c-axis contraction, by 
taking into account the c-axis compressibility^^ and the 
variation of J with pressure. Thus it is concluded that 
the lattice strains would affect J{x) and ctMg only to a 
minor extent. 

Now, we shall determine the variation of the spin- 



FIG. 4: Doping dependence of Tn{x) (top) and of 
M{x)/M{<S) = Bl{x,T 0)/Bl{0,T ^ 0) (bottom) in 
La2Cui_:^Mg^04 and in La2Cui-;i;Zn:^04 (B/,(0,r ^ 0) = 
427 Gauss). The dashed lines are guides to the eye, while the 
solid line in the bottom panel shows the behaviour expected 
for a diluted QJ1J2SL (see text). In the inset the ratio of 
the spin stiffness derived for Zn and Mg-doped compounds is 
reported. 



stiffness pf"'{x) due to Zn doping. Given the sim- 
ilar trend of the lattice parameters in Mg and Zn 
doped samples (Fig. 3), it is reasonable to use also 
for La2Cui_a;Zna;04 the values of J{x) derived for 
La2Cui_a;Mga;04. From Eq. |31 by taking the ratio of 
T/v (a;) for Zn and Mg-doped compounds one can write 



\x) 



rpZr 



{x),^f t(x) 



(x) 



2.3J{x)t{x){l-x)^ 
An 



where t{x) — 
M^"(x)/A//^'^f(a;) 



n^\x)/T^''' 



(4) 

By taking the experimental data for 
Tpf{x) and M{x), on the basis of EqlH one finds a system- 
atic reduction of pf"(a;) with respect to pf^^{x), which 
reaches about 16% in the x=0.12 sample (inset to Fig. 
4). Namely, there is a decrease of the spin-stiffness for 
Zn-doped samples which is much faster than the one 
expected for a dilution model. Since the reduction of 
pf"(a;) cannot be due to lattice strains, as mentioned 
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above, it is clear that the observed variation should be 
associated with a different effect of Zn and Mg orbitals 
on the superexchange coupling mechanisms. In particu- 
lar the marked reduction of the spin stiffness reported in 
Fig. 4 suggests that Zn doping gives rise to interactions 
that compete with J, namely that there is a frustration 
induced by the impurities, consistently with previous the- 
oretical predictionsji^ 

If now one considers the doping dependence of M{x) 
one finds a very similar trend both for Mg and Zn- 
doped samples (Fig. 4). Again, one can consider that 
La2Cui_:EMg^04 behaves as a diluted spin system with 
J{x) slightly increasing with doping due to the structural 
modifications. Also here it is instructive to observe that 
the lattice contraction associated with the application of 
an external pressure yields a small but non-negligible re- 
duction of Af(x).— Thus, it is tempting to associate the 
slight decrease of M{x) found in La2Cui_2:Mg2:04 with 
respect to the dilution model^^ to lattice effects. On 
the other hand, it is not clear why in La2Cui_j:Zna;04, 
where "(x) is reduced, M{x) does not show a more pro- 
nounced decrease. This implies that the more significant 
reduction of Tn{x) found in La2Cui_2;Zna;04 originates 
from the decrease of pf"(x) rather than from the one of 
M{x) (see Eq. 3). 

In order to understand why M{x) reduction is slightly 
lower in La2Cui_xZnj;04 one has to consider that if Zn 
is giving rise to a next n.n. coupling, the magnetic lat- 
tice around the impurity is no longer a square lattice^^ 
and the mapping of the microscopic Hamiltonian onto 
a 2DQIIAFSL Hamiltonian should no longer be valid. 
Indeed, in that case the local configuration around Zn 
would be more similar to the one found in the diluted 
S = 1/2 Ji — J2 model on a square lattice (QJ1J2SL),^'' 
where the more connected spin texture yields a less pro- 



nounced decrease of M{x). In fact, in Li2Vi_2;Tia;Si04)^ 
a prototype for the diluted QJ1J2SL model, the initial de- 
crease of M{x)/M{0) ~ 1 — 0.46x is lower than the one 
found for a diluted 2DQHAFSL (Fig. 4)." Still one has 
to consider the effect of "(a;) reduction on M{x). In the 
framework of the Ji — J2 modeP^ the observed reduction 
of pf " (x) would correspond to an effective increase of the 
ratio J2/J1 to~ 0.04 for x — 0.12, which should lead to a 
correspondent reduction of M{x)/M{0) by only 3.6%;^ 
i.e. to an additional decrease M{x)/M{0) ~ 1 — 0.3x. 
Thus, taking into account both the effect of dilution 
and the increase of frustration the observed decrease of 
M{x)/M{0) ~ 1 -0.75a; due to Zn doping (Fig. 4) seems 
reasonable. Qualitatively speaking, the microscopic con- 
figuration that one should find in La2Cui_j;Zna;04 is 
in between the one of a diluted 2DQIIAFSL and the 
diluted QJ1J2SL model and, accordingly, M{x)/M{0) 
curve should stay in between the trend expected for those 
two models (Fig.4), as it is observed. 

In conclusion, from the comparison of M{x) and T]\[{x) 
in Mg and in Zn doped La2Cu04, it was found that 
La2Cui_j;Mgj;04 can still be described in terms of a 
dilution model with minor corrections due to lattice 
strains. On the other hand, the marked reduction of 
the spin stiffness found in La2Cui_2;Zna;04 indicates that 
in this latter system competing next n.n. interactions 
may arise around the impurity and generate a frustrated 
magnetic lattice. Also the reduced dependence of M{x) 
found in the Zn-doped system suggests the presence of 
a spin texture more connected than the one of a diluted 
2DQHAFSL, which indicates a different effect of Zn and 
Mg orbitals in the superexchange coupling mechanisms. 

Useful discussions with A.L.Chernyshev and the tech- 
nical assistance by Sean Giblin during the measurements 
at ISIS are gratefully acknowledged. 
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